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Abstract 

In this article, an analysis is carried out to study the effects of aligned magnetic field, radiation 
absorption and viscous dissipation on the MHD unsteady convective heat and mass transfer flow of a 
viscous incompressible electrically conducting and heat absorbing fluid along a vertical porous plate 
embedded in a porous medium with variable temperature and concentration. Approximate solutions 
for velocity, temperature and concentration are obtained by solving the governing equations of the 
flow field using multi parameter perturbation technique. The expressions for skin friction, the rate of 
heat transfer and mass transfer from the plate to the fluid are derived in non-dimensional form. The 
effects of various flow parameters affecting the flow field are discussed. It is found that the velocity 
increases with an increasing thermal Grashof number or massGrashof number or permeability 
parameter or viscous dissipationwhereas it decreases for magnetic field parameter or heat 
sourceparameteror radiation absorption coefficient or angle <j) . The temperature increases for 
increasing viscous dissipation while it decreases for increasing Prandtl number or radiation 
absorption coefficient in the presence of heat source parameter. A growing heat source parameter or 
heat sink parameter retards the temperature. 
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1. Introduction 

Flow problems through porous media over flat surfaces are of great theoretical as well as 
practical interest in view of their applications in various fields such as aerodynamics, extraction 
of plastic sheets, cooling of infinite metallic plates in a cool bath, liquid film condensation 
process and in major fields of glass and polymer industries. The study of heat and mass transfer 
with magnetic field effect is of considerable importance in chemical and hydrometallurgical 
industries. Soundalgekar [1] studied the viscous dissipation effects on unsteady free convective 
flow past an infinite vertical porous plate with constant suction. The two dimensional unsteady 
free convective and mass transfer flow of an incompressible viscous dissipative and electrically 
conducting fluid past an infinite vertical porous plate was examined by Gregantopouloset al.[2]. 
Kinyanjui et al. [3] solved the problem of MHD free convective heat and mass transfer of a heat 
generating fluid pastan impulsively started infinite vertical porous plate with Hall current and 
radiation absorption by using a finite difference scheme. The effect of the viscous dissipation 
term along with temperature dependent heat source/ sink on momentum, heat and mass transfer 
in a visco-elastic fluid flow over an accelerating surface was studied by Sonth et al. [4]. 
Cookeyet al. [5] investigated the influence of viscous dissipation and radiation on unsteady 
MHD free convective flow past an infinite heated vertical plate in a porous medium with time 
dependent suction. Aissaand Mohammadein [6] have analyzed the effects of the magnetic 
parameter, Joule heating, viscous dissipation and heat generation on the MHD micropolar fluids 
that passed through a stretching sheet. Salem [7] investigated the coupled heat and mass transfer 
in Darcy-Forchheimer mixed convection from a vertical flat plat embedded in a fluid- saturated 
porous medium under the effects of radiation and viscous dissipation. Zueco [8] used network 
simulation method (NSM) to study the effects of viscous dissipation and radiation on unsteady 
MHD free convective flow past a vertical porous plate. Prasad and Reddy [9] investigated 
radiation and mass transfer effects on an unsteady MHD free convective flow past a semi-infinite 
vertical permeable moving plate with viscous dissipation. The effects of viscous and joules 
dissipation on MHD flow, heat and mass transfer past a stretching porous surface embedded in 
porous medium were studied by Anjali Devi and Ganga [10]. HemantPoonia and Chaudary[ll] 
haveanalyzed the heat and mass transfer flow with viscous dissipationon an unsteady mixed 
convective flow along a vertical plate embedded in porous mediumwith suction. The effects of 
thermal radiation and variable viscosity on the unsteady hydromagnetic flow of an electrically 



A Monthly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., ItMJiWBtffi fi as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Management, IT and Engineering 

http://www.ijmra.us 



May 

2013 




Volume 3, Issue 5 



ISSN: 2249-0558 



conducting fluid over a porous vertical plate in the presence of viscous dissipation and time- 
dependent- suction has been presented by Mahmoud[12]. Ahmed and Batin [13] presented an 
analytical model for MHD mixed convective radiating fluid with viscous 
dissipation.BalaSidduluMalga and NaikotiKishan [14] have studied the effects of viscous 
dissipation on unsteady free convection and mass transfer boundary layer flow past an 
accelerated infinite vertical porous flat plate with suction. The effect of magnetic field on 
unsteady free convective flow of a viscous incompressible electrically conducting fluid past an 
infinite vertical porous flat plate embedded in a porous medium in the presence of constant 
suction and heat sink has been studied by Das et al. [15]. 

In this article an attempt is made to study the effects of aligned magnetic field, radiation 
absorption and viscous dissipation on the unsteady convective heat and mass transfer flow along 
a vertical porous flat surface through a porous medium with heat source/sink. 



2. Mathematical Formulation 

The two dimensional unsteady free convective flow of a laminar viscous incompressible 
electrically conducting and heat (radiation) absorbing fluid past an infinite vertical porous plate 
embedded in an uniform porous medium in the presence ofheat source or sink with constant 
suction under the action of aligned magnetic field strength Bo has been considered, x -axis is 
taken vertically upward direction along the plate and y -axis normal to it as shown in figure 1. 
In order to derive the fundamental equations we assume that (i) the flow variables are functions 
of y and t only, since the plate is infinite in extent (ii) p the density of the fluidto be constant (hi) 

the magnetic Reynolds number is small so that the induced magnetic field can be neglected (iv) 
the Hall effect, electrical effect and polarization effect are neglected (v) theJoule's dissipation 
term in the energy equation is neglected(vi) due to the application of suction at the surface, the 
fluid particles at the edge of the boundary layer will have a tendency to get displaced towards the 
plate surface, therefore v —> — v at y — > oo and this phenomenon is clearly supported by the 

equation of continuity.By usingBoussinesq's approximation, the governing equations of the flow 
field are given by 



A Monthly Double-Blind Peer Reviewed Refereed Open Access International e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., [tMJiWBtffi fi as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Management, IT and Engineering 

http://www.ijmra.us 



392 



May 

2013 




Volume 3, Issue 5 



ISSN: 2249-0558 




Fig 1 The flow configuration and co-ordinate system. 



Equation of continuity 



dy 



= 0^>v =-v 



(1) 



Equation of Momentum 



du du d 2 u „.< . . . . aBlSin 2 d , v ■ 

— + v — = V—^+gj3(T -TJ + g/3 (C -CJ- °- — u - — u 

at ay dy p k 
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Equation of Energy 
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dt' dy~ ^ 

The initial and boundary conditions are 



u =0,v =-v , r' =7 w + f (T w -rj)e im * ,C = C w + e (C H -C m )e im * at y = 

w — > , r — > , C — > as y — > oo (5) 



Introducing the following non-dimensional quantities 



y =- 



,f =- 



t_V^ 

4v 



Avco u' u n T -T 

-,u = —,v = —,& = —, " 

v P T w -T„ 



,C = 



c -c 
c, -c' 



M 



p >% 



Pr = — £c 



(6) 



in equations (2), (3) and (4) under the boundary condition (5), we get 



1 du du d 2 



4 dt dy dy 

1 dd dd _ 1 d 2 6 

4 dt dy~ Pr dy 2 

iac_ i a 2 c 

4 3? ~5c 3/ 



+ Gr <9+Gm C- 



' 2 ^ 
MSin 2 0+- 



+ Ec 



4 4 



(7) 

(8) 
(9) 
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The corresponding boundary conditions are 



u=O,0=l + £e io * ,C=l + ee im at y = 
u — > 0, — y 0, C — > asy^>°° 



(10) 



Here gis the acceleration due to gravity, p is the density, o is the electrical conductivity, [3 
is the coefficient of volumetric thermal expansion, J3* is the coefficient of volumetric mass 

expansion, v is a constant suction velocity, v is the coefficient of kinematic viscosity, co is the 

angular frequency, [1 is the coefficient of viscosity, Kis the thermal diffusivity, T is the 

temperature, !T w is the temperature at the plate, 7^ is the temperature at infinity, c p is the specific 

heat at constant pressure, Pris the Prandtl number, Sc is the Schmidt number, M is the magnetic 
field parameter, kis the permeability parameter, Gr is the Grashof number for heat transfer, Gm 
is the Grashof number for mass transfer, Qis the heat source/sink parameter, g, is the radiation 

absorption coefficient , <j) is an align angle and Ec is the viscous dissipation or Eckert number. 
3. Method of Solution 

In order to solve equations (7), (8) and (9) we assumes to be very smallandthe concentration, 
temperature, velocity of the flow field in the neighborhood of the plate as 



C(v,r) = C (v) + f e " w C 1 (y) 
0(y,t) = Q (y) + ee iia l (y) 
u(y,t) = u (y) + £e" a u x (y) 



(11) 
(12) 
(13) 



Substituting equations (11) to (13) in to equations (7) to (9) respectively and equating the 
harmonic and non-harmonic terms and neglecting the coefficients ofe 2 , we get 
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Zeroth order equations 



C = 



(14) 
(15) 



„ 2 1 

u + u -(M 5m — )w =-Gr0 o -Gm C Q 
k 



(16) 



First order equations 



c;-^c=o 



3 + Prf^ + 



Q-ico 



-VrOC 

Prfl = ^ITL _ 2PrEc u n u 



o "1 



J 



u x + u x -(M Sin 2 (f>+—+—)u x = -Grd x - GmC l 
k 4 

The corresponding boundary conditions are 



(17) 
(18) 
(19) 



y = : u =0, O = 1, C =1, u x = 0, d x =\, C x =\ 
y _> oo : Mq = 0, O = 0,C =0, Mj = 0, X =O,C X =O 



(20) 



Equations (14)-(19) are non- linear differential equations. Using multi parameter perturbation 
technique and choosing Ec«l, we assume 



#o - ^00 + Ec #01 



— ^00 ^01 



e x = e w +E C e lx 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
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Now using the equations (21) to (26) in to equations (14) to (19) and equating the coefficients of 
like powers of Ec, neglecting those of Ec because Eckert number Ec is very small for 
incompressible fluid flows, we get the following set of differential equations. 



Zeroth order equations 



C 00 =0 



c _Sciw c =Q 

Mo A Mo u 



(27) 
(28) 



Prg 

#oo + Pr #oo + #oo '- ^ 



pre, 



00 



Pr „ -Prg, 



#io + Pr #io + (2 - m — # 10 - »- io 



C„ 



M oo + M oo 



M 10 + M 10 ' 



M Sin 2 + 



1 



v 



u m =-Gr e m -GmC m 



1 

M Sin d> + — + — 
k 4 



m 10 = -Gr 6 W - Gm C w 



(29) 
(30) 
(31) 
(32) 



The corresponding boundary conditions are 



y-0 : u m - 0,6 00 -1,C 00 - l,u w - 0,6 lQ -1,C 10 -1 
J -> 00 : "oo = °» #oo = 0» C oo = 0» M io = 0, 0,o = 0, C 10 = 



(33) 



First order equations 

C , =0 (34) 
C„-^C n =0 (35) 



#o'i + Pr #oi + ^#o, = -^Qi - Pr i 



(36) 
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n + Px0 n +(Q-iQ)) — U = — C ~ 2 Pi u,ji 



4 11 4 11 



OCTMO 



M 01 + U 0l 



M Sin 2 + 



1 



w 01 = —Gr m - Gm C 01 



2,1 ifi^ 



M 5/n <Z> + — + 

* 4 



w n = —Gr0 u - GmC u 



(37) 
(38) 

(39) 



The corresponding boundary conditions are 



y = ;u m =O,0 Ql =0,C 01 =0,u u =O,0 U =0,C n =0 
y : m i = O,0 O1 =0,C 01 = 0,m u = O,0 n =0,C n =0 



(40) 



The ordinary differential equations (27) to (32) and (34) to (39) aresolved subject to the 
boundary conditions (33) and (40) respectively. Then substituting the solutions in to equations 
(21) to (26), we obtained the exact solutions for concentration, temperature and velocity as 
follows: 



1 0)t —tn,y 



C(y,t) = l+e e im e 
0(y, t) = [( -A 3 ) + Ec {A 1% e 1 '^ + A 19 e' 2 '^ + A 3Q e 



(41) 



-(m 3 +m 5 ) y 



fA 31 e-^)]+ 



£ e ia * [{A, e~" h - y - A l e~' v ) + Ec(A 2l e 



-(m 4 +m 5 ) v 



' + A 22 e 



-(m 3 +m 4 )y 



+ 



A 23 e (w)y + A 2 , e - ( " h - +m > )y + A 25 e - (w)y + A 26 e - (m,+ " h)y + A 21 e" hy ] 



(42) 



u(y,t) =[( A 10 e'^ + A 13 + A 14 e** ) + Ec ( A 40 e^ y + A 4] e^ y + A 42 e~ ( ^ ,y + A 43 e^ y + A 44 )] 
+ ee i6 *[(A 5 e'^ + \e-'"' y + A, e~ m * y ) + Ec(A 32 e <m ^ y + A 33 e <m}y + A 34 e {m ^ )y + 



(43) 
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Skin-friction coefficient r at the plate is given by 
du 



(44) 



Using equations (43) and (44), we obtained the skin - friction at the plate in non -dimensional 
form as follows 

t = [(-ra 3 A 1() - m 5 A u ) + Ec(-2m 5 A 40 - 2m 3 A 4l - (m 3 + m 5 )A 42 - m 3 A 43 - m 5 A^ )] + 

ee'^li—m^ —m 2 A 5 — m 4 A)) + Ec(—(m 4 +m 5 )A 32 — (m 3 + m 4 )A 33 — (m 2 +m 5 )A 34 — (45) 
(m 2 + m 3 ) A 35 — (m t + m 5 ) A 36 — (mj + m 3 )A 31 — m 2 A 3& — m 4 A 39 )] 



5. Nusselt Number 

The rate of heat transfer coefficient Nu at the plate is given by 

•V. (46) 

Using equations (42) and (46), we obtained Nusselt number in non -dimensional form as follows 
Nu = \—m 3 A 4 + Ec(—2m 5 A 2S — 2m 3 A 2g — (m 3 + m 5 )A 30 — m 3 A 3l )] + 

£ e'°* [(m l Aj — m 2 A 2 ) + Ec(—(m 4 + m 5 ) A 2l — (m 3 + m 4 )A 22 — (m 2 + m 5 )A 23 — (47) 

(m 2 + m 3 )A 24 — (n\ + m 5 )A li — {m^ + m 3 )A 26 — m 2 A 21 )] 



Nu = 



dy 



6. Sherwood Number 



The rate of mass transfer coefficient Sh at the plate is given by 
'dc 



Sh = 



dy 



(48) 



y=0 



Using equations (41) and (48), we obtained Sherwood number in non -dimensional form as 
follows 

Sh=-m l £e ia * (49) 
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I 1 ico Pr Pr O 
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38 



(m 2 + m 5 ) — (m 2 + m 5 ) — M 



2 1 -fi 5 y i« 2 

-GrA, 8 



4 m 2 -2m, -AT, 



, A 



-Gr Ajg 



L 41 



4m, — 2m, — M, 



A2 _ 



-GrA 



-31 



- Gr Ao 4 = 

m 3 — m 3 — Mj 



(m 3 + m 5 ) - (m 3 + m 5 )-M l 



•5 1 

5 A4 = — ( Ao Ai A2 A3 ) 



7. Results and Discussion 



In order to get an insight in the physical situation of the problem, the numerical values of 
the velocity, temperature, concentration, skin friction, Nusselt number and Sherwood number at 
the plate are obtained for different values of the physical parameters involved in the flow field 
and are analyzed with the help of figures from 2 to 13. The value of 5c is taken to be 0.66 which 
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corresponds to water- vapor and Pr is taken to be 0.71 which corresponds to air at25°C 
temperature with one atmospheric pressure, cot is taken to be k 12 and 0is taken to be ;r/6.The 
values of the other physical parameters are chosen arbitrarily. 

The concentration profile is plotted for different values of Schmidt number 5c in figure 2. 
It is observed that the effect of increasing values of Sc is to decrease the concentration.The 
variation of the temperature distribution for various values of Pr, Ec, Q and Q l are represented in 

figures 3 to 6. Figure 3 shows the effect of Prandtl number Pr on the temperature distribution. It 
is obvious that with the increase in the values of Pr, the temperature across the boundary layer 
decreases. Figure 4 illustrates the influence of the viscous dissipation Ec on the temperature 
profile in the boundary layer with respect to heat source parameter. It has been observed that as 
Ec increases, the temperature increases. From numerical calculations the same trend is noticed 
in the case of heat sink parameter (Q<0). 

Figure 5depicts the dimensionless temperature profiles for different values of heat source 
parameter (<2>0).It is noted that increasing the values of heat source parameter causes a 
reduction in the fluid temperature. From numerical calculations the same phenomenon is 
observed in the case of heat sink parameter (Q<0). Figure 6shows the effect of radiation 
absorption coefficient Q l on the temperature field due toheat source parameter. It is observed 

that the temperature decreases as Q l increases. From numerical calculations a reverse trend is 

noticed in the case of heat sink parameter. 

Figures 7 to 13 indicate the variations of the non-dimensional velocity for various flow 
parameters respectively. In figure 7, the velocity profile is plotted for various values of thermal 
Grashof number Gr. It is observed that the main stream velocity increases with an increase in the 
thermal Grashof number Gr.From numerical calculations, the same trend is noticed with the 
effect of mass Grashof number Gm. The influence of the magnetic field parameter M on velocity 
profile is predicted in figure 8. It shows that a growing magnetic field parameter retards the 
velocity of the flow field at all points. The effect of permeability parameter k is studied and the 
results are exhibited in figure 9. It is observed that the velocity increases with 
increasingpermeability parameter k. Figure 10 depicts the effect of an angle <f> on the velocity 
field. The magnitude of the velocity decreases with increase of angle <f> . 
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The variation of the velocity profile is shown in the figure 11 for varying values of 
viscous dissipation £cwith respect to heat source parameter. The velocity increases with an 
increase of Ec. Figure 12 is a graphical representation which depicts the velocity distribution for 
different values of heat source parameter. It is clear that, increasing heat source parameter retards 
the velocity of the flow field at all points. The influence of radiation absorption coefficient g, on 
the velocity profile is illustrated in figure 13 in the presence of heat source parameter. It is 
observed that the velocity decreases with an increase of g, . 




Sc=0 6D.0.B6 ,078,0.94 



II 2 2.5 3 3.5 4 4.5 5 
y *■ 




Fig. 2 Effect of Sc on concentration 
field when 00=5.0, £=0.2 



Fig. 3 Effect of Pr on temperature 

field when Gr=5, Gm=l, M-l, 
£c=0.002, K=\, Q=l, Qi=0.5, 
8=0.2, ©=1.0 
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Effect of Ec on temperature 
field when Gr=10, Gm=10, 
M=\, k=\, Q=l, <2i=0.5, e=0.2, 



Fig. 5 Effect of Q (for heat source 
parameter) on Temperature 
field when Gr=5, Gm=5, 
Ec=0.002, M=l, k=l, <2i=0.5, 
e=0.2, (0=1.0 





Fig. 6 Effect of gi on temperature 
field when Gr=5, 

Gm=5,Ec=0.002 M=l, k=5, 
2=1,6=0.2, co= 1.0 



Fig. 7 Effect of Gr on velocity field 
when Gm=0.1, £0=0.001, M=l, 
fc=0.1,e=l,Qi=0.01, e=0.1, 
co=1.0 
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Fig. 10 Effect of (j) on velocity field 
when Gr=5, Gm-0. 1 , 
£c=0.001,M=l, k=0.1,Q=l, 
g y =0.01,£=0.1, (0=1.0 



Fig. 1 1 Effect of Ec on velocity field 
when Gr=5, Gm=0.1, M=l, 
k=\,Q=\, Q^O.01,^0.1, 
co=1.0 
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Fig. 12 Effect of Q (for heat source 
parameter) on velocity field 
when G?=5, Gm=0.1, M=l, 
fc=O.l,02=O.Ol,e=O.l, co=l, 
£■6=0.001 



Fig. 13 Effect of Q i on velocity field 
when G?=5, Gm=0.1, M=l, 
£=0.1,2=1,6=0.1, (0=1, 
£c=0.001 



8. Conclusions 



Based on the results and discussions, the following conclusions have been arrived at. Increasing 
the Schmidt number induces reduction in the concentration andrises the rate of mass transfer and 
skin friction coefficient. The velocity and skin friction increase for the increase of Grashof 
number for heat and mass transfer or permeability parameter and decreases for the increase of 
magnetic field parameter or Prandtl number or angle (j) . With the increase of heat source 
parameter; the velocity, temperature and skin friction decrease while the rate of heat transfer 
increases. Also for the increase of heat sink parameter the velocity and temperature 
increasewhereas the skin friction and the rate of heat transfer decrease. The velocity, 
temperature, skin friction and the rate of heat transfer decrease with increase in radiation 
absorption coefficient in the case when Q>0. But the trend is just reversed in the case when 
<2<0.The temperature, velocity, skin friction and the rate of heat transferincrease with the 
increase in viscous dissipation Ec in bothheat source and heat sink parameters. 
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